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Abstract. Wc study, on the basis of the general entangled-plaquette variational 
ansatz, the ground-state properties of the spin-^ antiferromagnetic Heisenberg 
model on the triangular lattice. Our numerical estimates are in good agreement 
with available exact results and comparable, for large system sizes, to those 
computed via the best alternative numerical approaches, or by means of 
variational schemes based on specific (i.e., incorporating problem dependent 
terms) trial wave functions. The extrapolation to the thermodynamic limit 
of our results for lattices comprising up to Jf = 324 spins yields an upper 
bound of the ground-state energy per site (in units of the exchange coupling) 
of -0.5458(2) [-0.4074(1) for the XX model], while the estimated infinite-lattice 
order parameter is 0.3178(5) (i.e., approximately 64% of the classical value). 



PACS numbers: 02.70.Ss, 05.50+q 

1. Introduction 

In the last decades quantum Monte Carlo (QMC) pQ has emerged as one of the most 
powerful numerical tools to study the ground-state (GS) properties of quantum spin 
systems on a lattice. It has been successfully applied, for example, to determine 
quantities of crucial interest such as the energy and the order parameter of the two 
dimensional (2D) antiferromagnetic Heisenberg model (AHM) on the square lattice [5] 
and, more generally, can be considered the method of choice to investigate unfrustrated 
(bosonic) models in any dimension. 

The situation is more complicated when the geometry of the system is frustrated 
as in the triangular or kagome lattice. In this case QMC techniques based on imaginary 
time projection, which are "exact" for unfrustrated systems, suffer from a sign-problem 
and approximations or alternative techniques have to be employed. The density 
matrix renormalization group (DMRG) [3] method, although not affected by sign 
instability, yields very accurate results only for one dimensional (ID) spin systems, 
being instead hardly applicable to general 2D problems due to the unfavorable scaling 
of the computer resources needed with the system size [4]. 

Projected-entanglcd pair states (PEPS) [5], the natural extension to higher spatial 
dimensions of the matrix product states variational family, which constitutes the 
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foundation of DMRG, can efficiently approximate GS's of local Hamiltonians [5], 
and have been successfully employed to simulate frustrated and bipartite quantum 
spin systems [7] [5] GO QUI E]. Their applicability, however, is restricted by the 
disadvantageous scaling of the computational cost to 2D models with open boundary 
conditions. Different approaches, which have been proposed, and applied in 2D 
[12j [TjJ] HU [15] HH [17j HU [19] , suffer from similar limitations and seem hardly usable for 
large systems [possibly three dimensional (3D) ones] and periodic boundary conditions 
(PBC). The idea of combining DMRG and Monte Carlo (20] [21] has opened the 
opportunity to design new variational algorithms for Tensor-Network based ansatze 
[22l [23] . The applicability of one of them, based on string-bond states, has been 
demonstrated even in the case of 3D systems [22] . 

In a recent paper we have introduced a class of variational states, namely 
cntangled-plaquette states (EPS), which provides an accurate ansatz to simulate the 
GS properties of quantum spin models [23]. We have shown E3] that the EPS trial 
wave function ( WF) allows to describe efficiently the GS of lattice Hamiltonians taking 
advantage of a simple variational Monte Carlo algorithm. Extremely accurate results 
have been obtained by us [24] and other investigators (25] [26] for several frustrated 
(fcrmionic) and unfrustrated models. 

In this paper we study, as a paradigmatic example of frustrated system, the 
spin- 1 AHM on the triangular lattice. Previous works, carried out by means of both 
QMC approaches with a fixed-node approximation and variational calculations based 
on different classes of states, which explicitly incorporate antiferromagnetic ^ Neel 
order, found values of the energy per site extrapolated to the thermodynamic limit 

between 0.53 and 0.5470 (27] Ell O EH EU [32]. Evidence of long-range 

order has been generally found with important differences, however, in the estimate 
of the order parameter. Our work, is motivated, therefore, not only by the interesting 
physics which can emerge from the GS of a quantum antiferromagnet but also by 
methodological aspects: it is in fact important, in order to asses its applicability to 
generic quantum spin models, to benchmark the general EPS ansatz against other 
variational choices specifically designed to tackle a particular problem. 

We carried out computer simulations on systems of size up to N = 324 spins 
(previous applications of EPS on frustrated systems have been performed on much 
smaller system sizes, up to N = 64) [24]. We found for TV = 36, in the case of 
the XX model, an error on our estimated energy relative to the exact one [33] as 
small as 6 x 10~ 4 and an energy per site extrapolated to the thermodynamic limit 
of -0.4074(1). At the Heisenberg point we obtained E^ PS (N = 36) = -0.55420(4), 
which differs from the exact value by ~ 1%, and Eq PS (oo) = —0.5458(2), with a value 
for the order parameter of 0.3178(5). 

2. The EPS wave function 

In this section we describe the entangled-plaquette ansatz offering also some detail 
concerning the procedure adopted to minimize the energy and compute GS physical 
observables. Given a collection of N spins \ arranged on a lattice, a generic trial WF 
can be expressed as: 




(1) 



3 



where the VK(n)'s are complex weights, n) = \n\, 712, . . . , Hjv) and rii = ±1 V 
i = 1,2,...,N is the eigenvalue of <r?. The variational principle ensures that the 
energy expectation value on a trial state is a rigorous upper bound of the GS energy, 
therefore, such a quantity can be evaluated by minimizing the trial energy with respect 
to the weights. To achieve this task one has to specify the form of W(n), which 
corresponds to the choice of a particular variational ansatz. A desirable ansatz is one 
able to capture the GS physics of a system by means of few variational parameters, 
which can be efficiently optimized on a computer. The cntanglcd-plaquette ansatz is 
based on the following idea: 

i) Cover the systems with M plaquettes in such a way that I sites labelled by 
n i.p, n 2,p, ■ ■ ■ , ni.p = np belong to the Pth plaquette. 

ii) Express the weight of a global spin configuration as a product of variational 
coefficients C" 1 , C" 2 , . . . , C"/ f in biunivocal correspondence to the particular spin 
(e.g., along the z axis) configuration of the plaquettes. Hence: 

M 

(n|V>} = W(n)= J{cy. (2) 
p=i 

In order to gain additional insight into the EPS ansatz let us discuss few 
exemplificative cases: 

a) Each site is a plaquette: the wave function is of the mean field form (i.e., 
correlations are neglected). 

b) Plaquettes comprise more than one site and each site belongs to a single 
plaquette: with this choice it is possible to obtain reasonable estimates of the GS 
energy (up to corrections which scale with the plaquette boundaries) and short-range 
correlations (of the order of the plaquette size) . 

c) Overlapping (i.e., entangled) plaquettes are used: this ansatz clearly produces, 
accounting for the entangled nature of a quantum GS, estimates of the energy and 
long-range correlations much more accurate than those computable with a WF whose 
weights are the product of non overlapping plaquettes |34| 135] as in a) and b) . 

A variety of variational ansatze, as illustrated in [35], have a straigthforward 
representation in terms of EPS. This is, for example, the case of the WF introduced 
by Huse and Elser [27] for the AHM on the triangular lattice which contains up to 
three-spin correlations and a phase factor that depends on the system geometry. The 
accuracy of this variational choice can promptly be improved via a general (i.e., with 
no specific phase factor) EPS WF based on larger or differently shaped plaquettes, 
hence, including additional correlations. In our calculations we start from the trial 
state described in a) (i.e., N — M and 1 = 1) where the weights and the initial 
spin configuration in the computational basis are randomly selected; we estimate, via 
Monte Carlo sampling, the expectation value of the energy as well as its derivatives 
with respect to the Cp^'s; the value of these coefficients is then optimized by updating 
all of them of a small step against the gradient direction. Once convergence is 
reached the size of each plaquette is increased, naturally introducing correlations in 
the ansatz, and the procedure iterated. Finally, obtained the optimal GS WF, the 
mean value of observables other than the energy can be evaluated by means of the 
Monte Carlo method as well. An exhaustive discussion on the particular suitablility of 
EPS to efficiently perform the minimization procedure previously outlined is reported 
in |24j . It is important to mention that the entanglcd-plaquctte ansatz, as well as 
the proposed numerical protocol, is applicable to systems in any spatial dimension; 
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moreover our method, as any purely variational (i.e., not based on imaginary time 
projection techniques) one, is sign-problem free allowing for the characterization of 
both unfrtustrated and frustrated systems. 

3. Results 

In the following, we show results obtained with EPS for the 2D AHM on the triangular 
lattice. Computer simulations have been performed on several system sizes comprising 
up to N = 324 sites. Our numerical findings arc compared with exact diagonalization 
(ED) results available for N = 36 [33] and estimates obtained with alternative 
computational approaches or variational ansatze [27l [2Sj [Ml EH Ell E2] . The AHM 
Hamiltonian is: 

<i,j> <i,j> 

where J xy = J z = 1, S" is the spin-i operator in the a direction acting on site i 
and the summations run over nearest-neighbour sites. The calculations were carried 
on in the S^ OT = Sf = sector, being the total magnetization along z a 
good quantum number. Therefore, during the simulation, we update the system 
configuration proposing a spin- flip for couples of spins with opposite S z . As a first 
benchmark we estimate the GS energy of the XX model (i.e., J z = 0, J xy = 1) with 
PBC. Numerical results, shown in tabic Q] refer to different size I of the plaquettes. 
For I = 4 the relative error of the EPS energy of the N = 36 system with respect to 
the exact result is smaller than 2%, and decreases down to 6.5 x 10 -4 when I = 16. 
Our estimate, Eq PS (N = 36) = —0.41069(3), compares favorably to those obtained 
with alternative trial WF's [29l [32]. For example the trial state introduced in [27], 
usually employed as a starting point for fixed-node based QMC calculations, leads to 
a value of the energy per site of about —0.4039 [29] (i.e., approximately 1.7% higher 
than the exact result). 

The extrapolation to the infinite lattice limit of the energies in table [1] has 
been performed by assuming the functional scaling E^ PS (N) = E^ PS (oo) + 

the dependence of the GS energy per site on the system size, as well as the fitting 
functions for different values of I, are reported in figure [T] Extrapolated energy values 
range from -0.3987(2) (I = 4) to -0.4074(1) (I = 16) and are indistinguishable, taking 



Table 1. Variational GS energy per site (in units of J xy ) computed, using 
the entangled-plaqucttc ansatz, for the spin-i antiferromagnetic XX model on 
a triangular lattice of size TV with PBC. Error bars are in parenthesis; the number 
of spins per plaquette is reported in the first column. The exact result for TV = 36 
is also shown for comparison 1331 . 



1 


N=36 


N=144 


N=324 


4 


-0.40426(5) 


-0.39954(5) 


-0.39875(5) 


6 


-0.40626(4) 


-0.40197(4) 


-0.40120(4) 


9 


-0.40956(3) 


-0.40609(3) 


-0.40545(3) 


12 


-0.41006(3) 


-0.40685(3) 


-0.40630(3) 


16 


-0.41069(3) 


-0.40793(3) 


-0.40737(3) 


Exact 


-0.410957 
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Figure 1. Size dependence of the GS energy per site for the spin-i 
antifcrromagnctic XX model on the triangular lattice with PBC. Different symbols 
refer to different plaquette size I (i.e., stars: I = 4, boxes: £ = 6, diamonds: I = 9, 
triangles: I = 12, circles: I = 16. Error bars, not shown for clarity, are smaller 
than the symbol size. Each dashed line is a fit to numerical data with the same I 
(sec text). 



into account the statistical uncertainty, from the results obtained when N = 324. 
For this size, therefore, the thermodynamic limit appears fully recovered. The EPS 
extrapolated energy is slightly lower than that computed by extrapolating exact results 
up to N = 36 (i.e., Eq D (N = oo) = —0.4066), illustrating the importance of studying 
system sizes larger than those treatable via ED techniques. It has to be mentioned 
that the small deviation shown by our data in figure [T] from a linear scaling, being the 
same regardless the size of the plaquettes, should not be due to a loss in accuracy of 
the EPS WF as the system size increases; more plausibly, it can be ascribed to higher 
order corrections in the series expansion of the energy as a function of jr . Moreover 
the energy estimates (for a fixed lattice size) do not show evidence of convergence with 
the plaquette sizes that we have been able to utilize in this work. This is not surprising 
(in fact the EPS ansatz is formally exact when the plaquette size is as big as that of the 
whole system) and simply means that the exact result, even if the energy differences 
for I > 9 are quite small, is not yet reached. As a consequence one can expect, on 
further increasing the plaquette size, an improvement of the GS energy. It is however 
remarkable, taking into account the methodological aspect of our paper, that our best 
variational upper bounds, obtained by using plaquettes of 16 sites, already provide 
extremely accurate energy estimates. 

Next, we discuss our findings for the GS of the Hamiltonian ([3]) at the Heisenbcrg 
point (i.e., J xy = J z = 1). GS energies per site are reported in table [21 The EPS 
GS energy computed for N = 36 differs from the exact one [33] by ~ 1%, while, 
calculations based on ordered WF's yield typically higher values of the GS energy 

between 0.5367 [27] and 0.5515 [31]. When N increases, the general EPS 

WF remains competitive with respect to accurate choices specifically adapted to the 
triangular lattice geometry. Particularly, our energies for N = 144 and TV = 324 differ 
from the best (to our knowledge) variational estimates [31] by less than 0.25%. 



6 



Table 2. Variational GS energy per site (in units of J xy ) computed, by using the 
entangled-plaquette ansatz, for the spin-i AHM on a triangular lattice of size N 
with PBC. Error bars are in parenthesis; N plaqucttes comprising 16 sites have 
been employed. The exact result for N = 36 is also shown for comparison 1331 . 

N E^ PS Exact 

36 -0.55420(5) -0.560373 
144 -0.54700(8) 
324 -0.5459(1) 



As a result of a finite-size scaling analysis of the EPS energies per site carried 
out by extrapolating the numerical data in table [2] (see figure [2]) we find Eq PS (oo) = 
—0.5458(2), which is in excellent agreement with the most accurate results derived in 
previous works [29l [31] . 

GS estimates of the spin-spin correlation function CorrQ PS (rMAX,N) = (So • 
S rMAX ) e ps computed at the maximum distance (taking into account the PBC) on the 
lattice are reported in table [3] The magnetic order parameter of the system defined 
as m 2 = limis[- > . oo Corr (rMAX,N) is a quantity of central importance, being directly 
related to the presence of long-range order. Wether or not the spin-^ AHM displays 
long-range order has been a long lasting-question, however, the positive answer, on 
the basis of analytical [36l EI] and numerical [27l [28l HH [301 EU [33l [38l [39] results, 
is commonly accepted. Estimates of the order parameter expressed as a fraction of 
its classical value range from ~ 0.4 [29j [38] to ~ 0.72 [30]. Our computed spin-spin 
correlations for N = 144 and N = 324 are indistinguishable (within the statistical 
uncertainty). Hence we can consistently assume m EPS = (Corr EPS (ymax , N = 
324))i = 0.3178(5) (i.e., - 64% of the classical result). This value of the order 
parameter, approximately 2% lower than that obtained with the same numerical 



-0.5450 
-0.5470 

£ -0.5490 

m 

o 
W 

-0.5510 
-0.5530 



-0.5550 



-i 1 1 - 



-i 1 1— 



_i i i_ 



_i i i_ 



0.000 0.001 0.002 0.003 0.004 0.005 

N-3/2 



Figure 2. Size dependence of the GS energy per site of the spin-i AHM on the 
triangular lattice with PBC. Error bars, not shown for clarity, are smaller than 
the symbol size. The dashed line is a fit to numerical data. 
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Table 3. Variational GS spin-spin correlation (in units of J xy ), at the maximum 
distance (with PBC) on the lattice. Error bars are in parenthesis; N plaqucttes 
comprising 16 sites have been employed. 



N 


Corr$ PS {r MAX ,N) 


36 


0.1294(1) 


144 


0.1011(2) 


324 


0.1010(3) 



approach for the AHM on the square lattice [24], is consistent with the long-range 
ordered nature of the system found in previous numerical studies [29l [37l [38l EHl [30] . 

4. Conclusions 

The general entangled-plaquette wave function has been employed to study 
the ground-state properties of the spin-^ antiferromagnetic XX and Heisenberg 
Hamiltonian on the triangular lattice. With this ansatz and a simple variational 
Monte Carlo algorithm for the energy minimization, we found extremely accurate 
energy estimates for finite systems, as well as, by extrapolating, in the thermodynamic 
limit. The accuracy of our results compares well with that provided by the best 
alternative methods or variational ansatzc, which, when possible, (as in the case of the 
triangular lattice), are adjusted, to improve the estimates, by incorporating physical 
features of the particular system of interest. It is crucial to stress that the trial state 
utilized in this work, although not including any problem dependent term, provides 
a reasonable estimate even of the system order parameter; therefore, the variational 
class of entangled-plaquette states can be regarded as a reliable option to describe 
the GS of general frustrated (i.e., not treatable with standard quantum Monte Carlo 
methods) physical models on a lattice with no need of any a-priori knowledge. 
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